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ABSTRACT 



A program is developed to ascertain the response of a 
frictionless water column to internal waves. Using density 
strata as input, the program selects ten frequencies at 
equal intervals within a spectrum of internal waves bounded 
above by the maximum Vaisala frequency and below by the 
inertial frequency. Ray paths within the column are plotted 
for these frequencies. The first ten normal modes for each 
frequency are computed. The first three modes are plotted. 



2 



TABLE OP CONTENTS 



I. INTRODUCTION 7 

II. THEORY 11 

A. DEVELOPMENT OF THE INTERNAL WAVE EQUATION II 

B. APPLICATION OF RAY THEORY 14 

C. NORMAL MODE APPLICATION 17 

III. DEVELOPMENT OF THE MODEL 20 

A. MAIN PROGRAM 20 

B. SUBROUTINE INTERP 20 

C. SUBROUTINE RAY 20 

D. SUBROUTINE RUNK 21 

»Ti-r >Tr< rpiTT-i -r\T^r\r*r) ^ f on 

±V • J. J.J.NLT liilD J ' 

A. LINEAR INPUT 23 

B. THREE-LAYER INPUT 23 

V. APPLICATION OF THE PROGRAM TO OCEANOGRAPHIC DATA — 28 

VI. CONCLUSIONS 45 

APPENDIX - COMPUTER PROGRAM 46 

LIST OF REFERENCES 58 

INITIAL DISTRIBUTION LIST 59 

FORM DD 1473 60 



3 



LIST OF TABLES 



Table I. Comparison of Eigenvalues from Analytical 

and Numerical Solutions 24 

Table II. First Page of Computer Output for 

Deep Station 31 

Table III. Normal Modes Obtained for Deep Station 37 

Table IV. First Page of Computer Output for 

Shallow Station 38 

Table V. Normal Modes Obtained for Shallow Station 44 



4 



LIST OF FIGURES 



Figure 1. Variation of Vaisala Frequency with 

Temperature Profile 9 

Figure 2. Illustration of Angle Between Horizontal 

and Phase Velocity 16 

Figure 3. Relation Between Phase and Group Velocities - 18 

Figure 4. Comparison of Eigenfunctions for First three 

Modes of Frequency 0.00684 25 

Figure 5. Comparison of Vaisala Frequency Plots by 

Analytical and Numerical Methods 27 

Figure 6. Plot of Density Variation and Vaisala 

Frequency for the Deep Station 32 

Figure 7. Ray Paths for First Three Selected 

Frequencies of the Deep Station 33 

Figure 8. Ray Paths for Third through Sixth 

Frequencies of the Deep Sration 34 

Figure 9. Modal Structure for First Three Modes of 

First Frequency at Deep Station 35 

Figure 10. Density and Vaisala Frequency at Shallow 

Station 39 

Figure 11. Ray Plots for First Two Selected Frequencies 

at Shallow Station 40 

Figure 12. Ray Plots for Third through Ninth Frequency 

at Shallow Station 41 

Figure 13. Modal Structures for First Selected 

Frequency at Shallow Station 42 

Figure 14. Modal Structure for Third Selected Frequency 

at Shallow Station 43 



5 



ACKNO¥LEDGEMENT 



I sincerely thank Dr. Jerry Galt, friend as veil as 
advisor, who not only presented the initial problem, packaged 
neatly, but who also spent many extra hours helping me unwrap 
it. To my wife, Janice, I give special thanks for aiding me 
in drafting the thesis and, particularly, for weathering the 
storm of neglect while I concentrated on the problem. 



6 



I. INTRODUCTION 



Internal waves exist in all oceans, most bays and lakes; 
they vary widely in amplitude, period, and depth (La Fond, 
1962). Recently, advances have been made in interpreting the 
major features of the oceanic internal wave field (¥unsch and 
Dahlen, 1970). 

Internal waves can be observed as density perturbations 
through a stable, stratified water column. Since temperature 
correlates with density in the majority of the upper ocean 
areas, internal waves often appear as oscillations about the 
region of the thermocline, exhibiting greatest amplitude and 
pn as e speed w n e i' e i o i/ii e x uio c I x n e x s xn e m c s x 3. n t c n s c . —no 
oscillation initiated within the water column is described in 
terms of the Vaisala or stability freeuency. This frequency 
may be derived as follows (Zalkan, 1966). 

Consider the force acting on a water parcel raised 
vertically from its equilibrium position: 

/^P 2 

F = g[^- (P+?^Az)]dV-p ^^zdv ( 1 ) 

^ c 

where: g = gravitational acceleration, 

^ = density, 

£iZ = vertical displacement, 
dV = minute volume of the parcel, and 
c = sound velocity. The first term on the right hand 
side is the restoring force due to buoyancy. The second term 
is the contributing fox'ce of compressibility. The equation 
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for the induced vertical motion is: 



= ZdV 



( 2 ) 



The harmonic solution to equation (2) is the Vaisala frequency 
defined as 



N(z) - ,/^ 



(3) 



For oceanographic applications the water is effectively in- 
compressible and the second term under the radical is 
negligible. Thus, 



The Vaisala frequency is greatest in magnitude 
thermocline is most intense (see Figure l). 

The Vertical motion 'Initiated by the intern 
ing through the water column is governed by the 
equation (Fjeldstad, 1933): 



where the 



1 "np o a 

^ .. ^ j ^ 

differential 



where: V = vertical water velocity, 

k = internal wave number, 

(J^ = internal wave frequency, and 

f = inertial frequency. Equation (5) is a hyperbolic 
wave equation whose characteristics describe the direction 
of propagation of internal wave energy. With the appropriate 
boundary conditions that no vertical flow exists at the sur- 
face or at the bottom, the internal wave equation satisfies 
the criteria for a Sturmian system. The imposed boundary 
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TEMPERATURE — VAISALA FREQUENCY 




FIGURE 1. Variation of Vaisala Frequency 
with Temperature Profile 



9 



conditions are satisfied by the simultaneous propagation of 
two waves having the same horizontal wave number, one with 
its vertical component of phase velocity up, the other with 
its vertical component down. The stipulation that these 
waves cancel at the boundaries leads to the structure of 
normal modes. These modes constitute a discrete set of wave 
numbers that are the eigenvalues of the Sturmian system. 

This thesis, using as a basis FJeldstad's classic paper 
(Fjeldstad, 1933), d eveloped a numerical program that assessed 
the response by a frictionless, incompressible water column 
to internal waves in a frequency band bounded above by the 
maximum Vaisala frequency and below by the inertial frequency. 
Station density strata was used as input. Frequencies were 
selected within the band and their ray paths were drawn. 

Also, the governing differential equation was numerically 
integrated and the resulting endpoints were refined to deter- 
mine the first ten normal modes. The first three modes were 
plotted. 

The program was first tested using density strata with 
known output. The program then treated data taken from ocean 
stations off the coast of Monterey, California in mid-August, 
1970. 
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II. THEORY 



A. DEVELOPMENT OF THE INTERNAL ¥AVE EQUATION 

For a frictionless, incompressible water column the 
governing equations (following Feldstad, 1933) are: 
Momentum 



<)t 



1 ^ 



+ g + S ~ =0 



Continuity 



du , 



Densi uy invaria'^ i 

dt 









+ U V + V ^ = 0 



(6) 

(7) 

( 8 ) 

(9) 

( 10 ) 



The X-axis is directed in the direction of wave travel, 
the Y-axis is in the horizontal transverse plane of travel, 
and the Z-axis is vertical, positive upward. U, V, and V are 
the velocity components of the water in the X, Y, and Z axes, 
respectively. f is the coriolis parameter. P is the pressure, 
The depth is assumed to be uniform. Z = 0 is the bottom; Z = h 
is the undisturbed upper surface. For a steady state 
condition , 

dU dV dV 



dt “ dt “ dt 



= 0 . 



Thus , 



dP = - 



^gdZ. 



( 11 ) 



( 12 ) 
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P can be represented as 



P = Po'^S J pdz, (13) 

where Pq is the undisturbed pressure at the surface. Let 

P=(>q(Z) + P^(X,T,Z,t), (14) 

where is the undisturbed density and the perturbation. 
Therefore, pressure for the perturbation is: 



where 




(15) 



(16) 



Envoking a 
of density 
and squares 



general assumption that velocities and perturbations 
and pressure are small enough that their products 
may be neglected, equation (iO) riUiy be written: 



1 , 1 



V = 0 



^ 6 ^ 0 ^ ^ 



(17) 



0 



If, for the internal wave form, harmonic oscillations 
such as the following are considered: 

,P^ ^ cos (CTt - kx ) exp ( ay ) 

V,V sin (Cf^t-kx) exp (ay) 

where the constant a controls the transverse envelope and 

allows for the possibility of edge waves. The sign of a is 

such that the wave amplitude decays away from the boundary. 

The governing equations now become: 

P, 

(f' U - fV+ k^^= 0 (18) 



(T 



V + f U + ag^ = 0 
^0 



(19) 
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J 

I 



( 20 ) 






kU + aV + 4? = 0 

a/j 



( 21 ) 



Solving equations (l8) and (19) for U and V: 



U = 



k + f a 1 



CT2_f2^ 



0 



V = 



fk+ a ^1 



'0'2-f2 Qo* 

Solving for the vertical pressure gradient: 

^ q 



( 22 ) 



(23) 



dZ 






‘0 



V 



(24) 



where N = 



2 ^0 

~r^ as previously defined. When one inserts 

!o dZ ^ 



v22j and (2j) in v9) 



_ (O' ^ + f^) ^ 

2 2x dZ* 



When is eliminated from (24) and (25) and the transverse 
boundaries are disregarded (i.e., set a = 0), 

__MJ\a_W r 

A7 ' A'J I- ,o 



(25) 



d^V ^ 

,2 



dZ" 






(26) 



Equation (26) is the internal wave equation describing the 
response by a water column to a plane internal wave of known 
f requency , 0^ , The order of magnitude of the second term is 
small compared to the first and third. However, the numerical 
program can include the second term without any difficulty. 

It will be so retained. 
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B. APPLICATION OP RAT THEORY 



An analytical solution to equation (26) may be made by 
the transformation: 

V= (Z) exp ^ |). (27) 

¥hen the effect of the second term in eq. (26) is small, the 
analytical solution to the internal wave equation can, with 
a high degree of accuracy, be represented by'f* , the 
oscillatory portion. The oscillatory solution is: 






+ (- 



It 



-)(N^ _5J^)V= 0 



- f^ 

Solutions to equation (28) are of the form: 



( 28 ) 



=A exp i (k,X+k, ^ 5 - Y+^t) (29) 

^ ^ _ f^ 

where kj^ = horizontal wave number. 

Define© as the angle between the horizontal and the 
phase velocity. Then, 



c os 0 = 



= [ 



^2 _ ^2 1/2 



(k, + k, [ 



2 rN ^ '■^2 _ ^2 



] 



(T2_ 



■]) 



(30) 



Since 



o p 1/2 p p 1/2 

(^^-f^) =(N^-Q^^) 



(k,^ +k.^)^/^ 



where 



V =[- 






1 2 . , 2 
ki +k2 



-] 



n2_ 2 1/2 

^2 ^1 2 _ ^2 ^ 

|k| = (k 2 + k.2)^/^ 



(31) 



(31a) 



(32) 

(33) 
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Figure 2 illustrates the angle between the horizontal 
and the phase velocity. The phase velocity, broken into X 



and Z components, is; 

^1 

, 2 . , 2 " 



^1 r, 2.,2^, 2„2-|l/2 

2^1 2x3/2 '-1 ^ ^ ^ 



and 



^ + k2 ( k^ ^2 ^ 



^Z = F" ' " ' ' 2 372 [k^V + k 2f2]l/2 

^ (k^2 ^ j^^2)3/2 1 2 



For the group velocity. 



Cg =V, ^ 



where 



C 



S7 { ^ ^ \ 

Vk=<3k7 ^ J ) 

2,, 2 ^ ,, 2^2 1/2 



f SI r, 

"rr - l^L \ 



k "N"- + k^"*f 







The X-component of group velocity is 



2..2 . 2„2 1/2 



X k/N^+k^^f^ 
('^g)x=5k^[( + ^ 



^1^^2 r (N^ - f^) 1 

L“ o ^5 5 o . 1 /o J 



|kp (k^^N^ + 

The Z-component of group velocity is: 

2,, 2 ^ 2^2 1/2 



X kFN‘' + k^‘'f‘ 

'z = 5k7t 



(Cg)? = ^ — ) ] 



1 2 ^ , 2 
ki +k^ 



^1^^^2 r (N^ - f^) 



|kl 



[■ 



^ '■ (k^ V + k2"'f'^) 



2„2xl/2 



] 



Define© as the angle made with the horizontal by 
g 

group velocity vector. Then, 



(34) 

(35) 

(36) 

(37) 

/' O o \ 

(39) 

(40) 
the 
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FIGURE 2. Illustration of Angle between 
Horizontal and Phase Velocity 
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( 41 ) 






z 

X 






2 1/2 

2 ] 



The ray trace part of the program drew the ray paths as 
depicted in equation (41) for ten frequencies selected within 
the iner tio-gravity band. 

Determining the dot product of the phase and group 
velocities : 



C 0 C 

g 



k, (k,^N^ + k.^f^)^/^ k,k, 

( ^ V ^ )( ^ ^ 



|k| 



(N^ - f^) ^ 

O O O 1 /o / 



|kp (k^^N^ + k^^f^)^/^ 



k (k, V + k^^f^)^/^ k,^k„(N^ - f^) 

■ ‘ lip ’ ‘ |Ppo/v7Ip?p72> 



= 0 



Thus the phase velocity is perpendicular to the group 
velocity. Also, 






k2(k2^ - k^^)f^ 



(k^^N^ + k2^f^)^'^^ 



(k,2^k77^ 



(42) 



(43) 



The phase velocity and group velocity then have a relationship 
as seen in Figure 3. 



C. NORMAL MODE APPLICATION 

With the use of suitable boundary conditions, equation 
(26) may be solved to obtain modal structures. The applicable 
conditions are: 

V=0 at Z = 0 (44) 
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iriiiiiii! 1 
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z 

A 




Vz 



Figure 3. Relation between Phase and Group Velocities 



18 



and 



dZ 



(44a) 



(T 




gV = 0 



at 



Z = h 



To a high degree of approximation the vertical velocity shear 
is negligible at the surface and the second condition becomes: 

V=0 at Z = h (45) 

Equations (26), (44), and (45) compose a Sturmian system. 

V, the vertical velocity, is the eigenfunction. For a given 
frequency, discrete wave numbers or eigenvalues are determined. 
The program utilized a fourth-order Runge-Kutta numerical 
method to obtain an end value after integration of the 
internal wave equation over the water column. The endpoint 
was refined using a Newton-Raphson iteration procedure. 
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III. DEVELOPMENT OF THE MODEL 



A. MAIN PROGRAM 

The main program prints the station density stratifica- 
tion. The program then calls upon the three sub-programs to 
interpolate density between station depths, to plot ray paths 
for selected internal wave frequencies, and to print out and 
plot modal structures for these frequencies. 

B. SUBROUTINE INTERP 

Subroutine INTERP interpolates the density stratification 
for 400 incremental depths through the water column. The 
interpolation is accomplished by the use of a third-degree 
polynomial that matches points and slopes at each depth where 
data was obtained from a Nansen bottle. To guarantee density 
stability, the smaller of the slope above or below each point 
is chosen as the slope at that point. From the interpolation 
polynomial, INTERP computes the Vaisala frequency for each 
density value. 

C. SUBROUTINE RAY 

Subroutine RAY plots the interpolated density and the 
Vaisala frequency as a function of depth. RAY then determines 
the inert io-gravity frequencies at rates of l/lOO, l/lO, 2 /l 0 , 
3/10, 4/10, 5/10, 6/10, 7/10, 8/10, and 9/10 of the frequency 
band. The program then plots the characteristic paths over 
the depth of the column for each of these frequencies. 
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Additionally, subroutine RAY determines the percentage of 



the water column depth where the selected frequency exceeds 
the Vaisala frequency. This percentage represents that por- 
tion of the column where the Sturmian solution to the internal 
wave equation is exponential. A greater percentage implies 
increasing difficulty in applying boundary conditions during 
the numerical integration. The percentage figure for each 
frequency is listed beneath the ray path plot. 

D. SUBROUTINE RUNK 

Subroutine RUNK utilizes a fourth-order Runge-Kutta 
numerical method to integrate the internal wave equation over 
the water column. The endpoint of the integration is refined 
by a ^n-Raph son iteration method until the endpoint is 

1/200 of the maximum value of the eigenfunction, V, determined 
during the integration scheme. For each selected frequency, 
the first approximation for the initial eigenvalue is 



The approximation is refined until the wave number correspond- 
ing to the first mode is determined. The approximation for 
the second mode is twice the value of . The initial esti- 
mate of the remaining eight eigenvalues is the previous 
eigenvalue plus the difference between the previous two 
eigenvalues . 

After the first three eigenvalues are determined, the 
program plots their modal structure. The remaining seven 
wave numbers are then listed. 






( 45 ) 



K 
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As the selected frequencies within the inertio-gravity 
band increase in magnitude, the boundary constraints, though 
analytica iiy well-defined, weaken exponentially. The higher 
frequencies are greater than the incremental Vaisala frequency 
for a greater percentage of column depth and it becomes 
increasingly difficult for the numerical program to find 
eigenvalues. RUNK terminates its search for the next eigen- 
values by (l) stopping after 20 refinements by the Newton- 
Raphson iteration process, (2) restricting the magnitude of 
the eigenfunction to 10 cm/sec, and (3) stopping if the 
solution converges to a mode other than the one searched for. 

Also, as the frequency within the band increases, V 
becomes more sensitive to changes in k when seeking higher 
modes. For situat-Lons where a minute change in k causes 
abrupt responses in end values of V, the eigenvalue may be 
only refined to a certain degree. These eigenvalues are 
indicated as approximations. 
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IV. TESTING THE PROGRAM 



A. LINEAR INPUT 



To check numerical formulation, coding, and accuracy, two 
test cases were used whose analytic solutions were possible. 
The first test was a linear variation of depth from 1.023 
grams/cc at the surface to 1.033 grams/cc at the bottom depth 
of 1000 meters. The Vaisala frequency for this input was 
constant and equal to: 



_ [ITK 



- 9. . 8 _ MZse _ c_ ^ 5=9.77X10 ^sec ^ 



(46) 



'1.025 614/cc 

The plot by the program of tiie deiiiit^ matches that of t-lie 
linear variation. The Vaisala frequency determined numeri- 
cally matches that determined analytically. Table I compares 
the analytical and numerical results of the linear variation 
of density for frequency 0.00684 cycles/min. Figure 4 
compares the first three modal structures. Analytically, the 
normal modes are 






K = ( o o 
m j^2 _ ^2 



The eigenfunction is 



/ N^ Z \ . / mlT „ N 

exp (— j) sin Z) 



(47) 



(48) 



B . THREE-LAYER INPUT 

The second input was a density strata that was constant 
at 1.023 grams/cc from the surface to 300 meters, varied 
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Table I. A comparison of the first ten normal modes 



k. 



k. 



^10 



for frequency 0. 


00684 cycles/min: 




Analytical 


Numerical 


% Erro: 


3.613 X 10"^ 


3.629 X 10“^ 


0.45 


7.225 X 10"^ 


7.259 X 10"^ 


0.47 


1.083 X 10""^ 


1.088 X 10“^ 


0.27 


1.445 X 10“"^ 


1.451 X lO”"^ 


0.42 


1.806 X lO""^ 


1.814 X 10""^ 


0.44 


2.167 X 10""^ 


2. 177 X 10"^ 


0.46 


2.529 X lO""^ 


2.540 X lO""^ 


0.43 


2.890 X 10”"^ 


2.903 X 10“"^ 


0.45 


3.251 X lO-"^ 


3.266 X lO""^ 


0.46 


3.613 X 10“"^ 


3.629 X 10“"^ 


0.44 
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linearly from 300 meters to a value of 1.027 grams/cc at 700 

meters, and then remained invariant from 700 meters to the 

bottom at a depth of 1000 meters. The corresponding value 

of the Vaisala frequency was then zero in the upper and lower 

-3 

isopycnal layers and equal to 9.8x10 in the middle region. 

The density plot by the program matched the analytic 
solution to the third decimal place. Figure 5 illustrates 
the comparison between the analytic solution and the computer 
solution for the Vaisala frequency, and illustrates how the 
interpolation scheme handles discontinuities in the density 
slope . 
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-1 




FIGURE 5. 



Comparison of Vaisala Frequency Plots 
by Analytical and Numerical Methods 
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V. APPLICATION OF THE PROGRAM TO OCEANOGRAPHIC DATA 



The program utilizes density strata from Nansen casts. 

It can process several stations. IDATA, the number of sta- 
tions to be studied, is the first card to be typed. It 
follows the alphanumer ics that are read in and serves in the 
plotting portions of the program. IDATA is read in as an 
integer, right-oriented to column ten of the IBM card. The 
second card to be typed up consists of three figures. The 
first is M, the number of bottles at the station, read in as 
an integer, right-oriented to column 10. The second is ZMAX, 
the depth at which the last bottle was taken, read in the 



format FlO.j. right-orieirlod lo coluipi- ?0. 



mi. „ 4-U A -P A 

i. I I \ . oliJxVu - 






is RLAT , the latitude of the station, that is read in using 
format PlO.l right-oriented to column 30. 

The next cards to be typed are those representing the 
density strata. These are read in as four pairs of numbers 
for each card. Each number has a field width of ten and is 
right-oriented. The first number in each pair represents 
station depth and has a format of FlO.l. The second of the 
pair has a format F10.5 and represents the station density. 
These four pairs of numbers, representing 80 columns per card, 
are entered until all the station data has been read in. 

The above represent all the cards needed to be typed for 
each station. From experience gained using the IBM 360 
computer to process ocean stations taken in mid-August, 1970, 
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off Monterey, California, it appears that two minutes are 
required to process shallower stations (around 200 meters) 
and four minutes to process deeper stations (around 2000 
meters ) . 

The actual printed output for each station appears as 
follows. On the first page of printed output appears: 

(l) the number of station depths, (2) the lowest bottle depth, 
(3) a printout of the station depths and densities in columnar 
form, ( 4 ) the inertial frequency computed for the station, 

( 5 ) the maximum Vaisala frequency computed for the station, 
and (6) the t en frequencies within the inert io-gravity band 
selected for consideration. 

The second page of printed outj)ut plots the Vaisala 
frequency and density versus depth. 

The next ten pages plot the ray path for each frequency 
versus depth. At the bottom of each ray plot, the amount of 
time that the frequency is greater than the Vaisala frequency 
is indicated. 

The remaining twenty pages of output display information 
on normal modes. Corresponding to each of the ten selected 
frequencies are two pages. The first page lists the first 
three normal modes and plots these modes versus depth. The 
second page lists the remaining seven modes. 

Two stations, one shallow and one deep, were chosen to 
illustrate the output from the program. The results for both 
stations are summarized by the following tables and figures. 
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Table II lists the information for the deep oceanographic 
station displayed on the first page of computer output. This 
station is representative of an oceanic internal wave regime. 
Figure 6 shows the plot of the density variation and the 
Vaisala frequency for the deep station. The Vaisala plot 
follows the pattern of the density profile. 

Figure 7 depicts the ray paths for the first three select- 
ed frequencies. The first frequency selected, 0.01371 
cycles/min, never exceeds the incremental Vaisala frequency. 
Thus, theoretically, its solution is always oscillatory and 
it travels horizontally nearly 30 km. The second frequency 
selected, 0.12969 cycles/min, exceeds the Vaisala frequency 
in the region of the water column below a depth of 970 meters. 
Thus it exponentially ’decays' in this region and its ray 
path over this portion of the water column is vertical. 

Figures 7 and 8 illustrate that, with increasing frequency, 
the horizontal distance covered by the ray path decreases 
and a greater portion is in the exponential region. The ray 
paths for the last four selected frequencies were not shown 
since their ray is nearly vertical from top to bottom. 

Figure 9 plots the first three modes for frequency 0.01371 
cycles/min. This was the only frequency for which the first 
three mornal modes could be positively identified. Besides 
gaining a better understanding of the actual appearance of 
the modal structures. Figure 9 suggests optimum placement of 
a thermocline follower for mode detection. For example, the 
follower would best detect mode three at a depth of about 
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TABLE II. Data display for deep station 



Incremental depth Incremental density 



0.0 

10.0 

20.0 

30.0 

50.0 

75.0 

100.0 

150.0 

200.0 

300.0 

400.0 

500.0 
600.0 
800.0 

1000.0 

1200.0 

1500.0 

2000.0 

2500.0 

3000.0 



1.02498 
1.02505 
1.02522 
1.02565 
1.02588 
1.02611 
1.02624 
1.02646 
1.02656 
1.02679 
1.02681 
1.02702 
1.02711 
1.02730 
1.02743 
1.02750 
1.02759 
1.02770 
] .02774 
1.02776 



Frequency niimbor 



Cycles/ min 



inertial 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

maximum Vaisala 



0.00082 

0.01371 

0.12969 

0.25855 

0.38741 

0.51628 

0.64514 

0.77401 

0.90287 

1.03173 

1.16059 

1.28946 
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FIGURE 6. Density variation and Vaisala frequency for deep station 
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FIGURE 7 . Ray patlis for deep station 
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FIGURE 8. Ray paths for third through sixth frequency at deep station. 
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FIGURE 9. Modal structures for first three modes of first frequency 



2000 meters. Table III lists the various modes determined 
for the selected frequencies for the deep oceanographic 
station . 

Table IV displays the first page of printed output for 
the shallow station representative of the regime in a large 
bay. Figure 10 illustrates the density variation and the 
Vaisala frequency as a function of depth. As in the deep 
station, the first selected frequency of the shallow station 
has an oscillatory solution throughout most of the water 
column. Its ray path travels horizontally nearly 2500 meters 
Its exponential region may be seen as a vertical path at a 
distance of 1955 meters. Figure 11 shows the effect of an 
increase in frequency without a considerable increase in the 
effect of exponential decay. The second frequency selected, 
0.15818 cycles/min., has an exponential solution over the 
water column in nearly the same region as the first. However 
frequency 2 only travels a horizontal distance of 230 meters. 
Figure 12 depicts the decrease in ray path as the frequency 
increases . 

Figure 13 depicts the first three normal modes for 
frequency 0.01656 cycles/min. Figure 14 illustrates the 
first three normal modes for frequency 0.31554 cycles/min. 
Figure 14 illustrates well the manner in which the modal 
structure deteriorates with increased frequency. Table V 
lists the modes determined for selected frequencies at the 
shallow station. 



36 



Table III. Normal modes for selected frequencies at 



the deep station. 



0. 01371. cycles/min. 

1. 0.00011736 

2. 0.00023344 

3. 0.00032185 

4. 0.00043948 

5. 0.00055503 

6. 0.00063828 

7. 0.00079211 

8. 0.00085791 

9. 0.00096286 

10. 0.00011641 



0. 12969. cycles/min. 

1. 0.0018469 

2. 0.0034051 
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Table IV. Data display for shallow station 



Incremental depth Incremental density 



0.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

40.0 

50.0 
60.0 
75.0 

100.0 



1.02535 

1.02537 

1.02543 

1.02551 

1.02570 

1.02580 

1.02606 

1.02610 

1.02614 

1.02614 

1.02616 

1.02619 



Frequency number 



Cyc les/min 



inertial 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

maximum Vaisala 



0.00083 

O AT At; A 

^ ♦ V-/ j v> ^ X./ 
0. 15818 
0.31554 
0.47289 
0.63024 
0.78760 
0.94495 
1.10231 
1.25966 
1.41702 
1.57437 
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FIGURE 10. Density profile and Vaisala frequency at the shallow station 
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FIGURE 11. Ray plots for first two selected frequencies at the shallow station. 
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FIGURE 12. Ray plots for third through ninth frequencies at the shallow station 



MODE 1 ^ODE 2 MODE 




42 



FIGURE 13. Modal structure for first selected frequency at shallow station 
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FIGURE 14. Modal structures for third selected frequency at shallow station 



Table V. Normal modes for selected frequencies at the 



shallow station. 



0. 


01656 cycles/min. 


0. 


15818 cycles/min. 


1. 


0.00087701 


1. 


0.0086086 


2. 


0.0025324 


2. 


0.025461 


3. 


0.0035674 


3. 


0.041537 


4. 


0.0050806 


4. 


0.051653 


5. 


0.0055217 


5. 


0.069448 




31554 cycles/min. 


6. 


0.074546 


0. 


7. 


0.092209 






8. 


0.10860 


1. 


0.018642 


9. 


0.12385 


2. 


0.056013 


10. 


0.13605 


3. 


0.10871 






4. 


0.14614 






5. 


0.17507 






6 . 


0.20301 






7. 


0.23839 






8. 


0.28283 






G 

y • 


0.33293 






10. 


0.37399 






0. 


47289 cycles/min. 


0. 


63024 cycles/min. 


1. 


0.032306 


1. 


0.053539 


2. 


0.094660 


2. 


0.14805 


3. 


0.19395 




78760 cycles/min. 


4. 


0.26040 


0. 


5. 


0.33316 






6 . 


0.40811 


1. 


0.092463 


7. 


0.48565 


2. 


0.23132 


8. 


0.55884 






9. 


0.63999 






10. 


0.72820 
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VI. CONCLUSIONS 



The numerical program developed in this thesis has derived 
a technique that links the theoretical application of internal 
waves to real data. Real solutions are determined for actual 
density profiles. The program can yield a useful description 
of the internal wave regime within a given area. Prom 
representative samples of stations investigated, four points 
are apparent; 

(1) , that with greater depth the effect of the boundary 
conditions weakens; 

(2) , that for low frequencies that exceed the Vaisala 



frequency over only a small percentage of Ihc water 
the ray path approaches the horizontal, while for h 
frequencies the ray path approaches the vertical; 

(3) , that the plot of the lowest three modes demons 
the manner in which, with increasing frequency, the 
of the boundary conditions weakens as the oscillate 
regions give way to exponential "decay" regions; 

(4) , that for the stations investigated by this pro 
only those frequencies within the lower quarter of th 
inertio-gravity frequency band could be solved for no 
modes with any degree of success. 



c o i umii 
igher 

trates 

effect 

ry 

gram, 

e 

rmal 



} 



45 



APPENDIX 





-If •- 






* 


1 


If — 








e Q 


If 0 








UJ LU 


If 0 






* 


3 of" 


If stiU 


0 






auj >-< 


If — 'X 


2 




if- 


UJXQO*-' 


If xa: 


— • 


2 


■Jt* 


QCH-UJQIO 


if < X 


1 — 


X 




LL h-LUa 


a- 001— 


1— 


X 


il- 


202:00 


If 


0 


X 




LLh-iUJUJoO 


•if <10 


0 




If 


OX-J < 


If >12 


X 


00 


If 


H-UJO 


If -f- 




< 


If 


Q»-hoo2oO 


If H- 


UJ 


X 


If 


23: ^LU 


If o^iu 


X 




If 


< XH-O 


If 02 


H- 


<I 


If 


CQ>0<0 


If 00 0 




h- 


If 


o<05: 


>f <J 


2 


< 


If 


<2LLI< 


If >x 


HH 


Q 


If 


LU CLO 


If 


X 




If 


23c^a< 


If x< 


0 X 


X 


If 


*-HOOa:x 


lf< ♦• 


CU uC 


0 


If 


xlullclq: 


If 2< 


00 X 


l-H 


If 


Ko: 0 


•k MO 


0 H- 


X 


If 


—•IX. «»X2 


If --x 


»« 




If 


s >0 


If --2 


X u 




If 


<0 2 


If 00 


cC 2 


SI 


If 


0002x1x1 


If 0 


< H- 


a 


Jf 


LU<XIKf- 


-iC- 




X 


If 


*~*003< 


If —*0 


2 X 


X 


If 


0»— OCLf- 


If 00 


' 2 




K- 


2<LU CO 


•'f XnT 


Q 


00 


If 


LU>o:>-a: 


If 


C 


2 


•if 


— ^ *X ^ ’ 


I< ^x 


CL 


0 


)f 


02 QTX 


•>f 


X 2 




If 


X130 


If ox 


X < 


h- 


If 


ci:2<xiLu 


If rOLU 


•> 


< 


If 


X 1“^ X X 


-u- 


X UJ 


H- 


If 


XHf-H- 


If- - 


CD 


00 


If 


2<a: 


If 1 — • 


O 




If 


Uj2 U.J0OOO 


If 00 0 


0 


X 


•jf 


21 — LU 


K- O-H 


h- a: 


0 


If 


LU-<02 


If X'— ' 


< 






oox -J— • 


If o:o 


00 \~ 


X 


If 


V— J— LUCL2 


If --UJ 


0 00 


Ul 


If 


0 X a: 


ij. —.X 


l-H 


CD 


i; 


X>-h~ 2iU 


K*ox 


OC 1 


2 


If 


OCQ <f- 0 


If •- 


X OrH 


0 


If 


XI >-XlU> 




2 Q< 


2 


If 


OOUICQOQO 


If M< 


0 *"4- 




If 


>02 


If *-h- 


2 ^ 


X < 


If 


2 0 3:tX2Lu 


if OOO 


< 2v4- 


X H- 


If 


<o:)Oa<3 


1C 


X 


H~ <4; 


If 


a:<-j o:o 


If H-X 


X xc\j 


0 


If 


0 LU LUOXI 




0 X 


0 — 


If 


oocnxoQ^ 


If OQ 


<< f~4 


f- — 0 


If 


0:X) H-Ct^x 


ifo: 




2 »-4 


If 


Q-QO a. 


If *-H- 


0 ^ 


LlJ f— 4 H-< 


If 


2^>- X 


If 20 


2 vO-- 


00 


If 


23<OLJO 


If 0 


KH 0 •* 


X in 


If 


•-•a 2X< 


If 2 


3:ooLnh- 


X 


If 


<0Q2xf-X 


■K- Oid: 


OX!— '< 


X 0 < 


1C' 


2 ^0 


If 22 


xx‘:a2 


xOQ2 


If 


oOOOOX 


If 2< 


xo<oc 


xxcx 


If- 


XX— J 1x00 J— 


If 00 


OCJXO 


XLUO 


If 


X^-^O 


If OCU 


X U J X X 


<xxx 


If 


f-OOXCS 


If If 


0 


H-X 


If 




l< 


xo 0 


<h- 0 


If 




If 


XX vO 


Q<t '-i 


If If If If If If If 1C- If If If 1^- 


h-CL 


•—•0 


ooooooououuo ooo ooo 







00 


Q 








< oo> 


Q 






X 




X 






K- 




X 






3 


<x<:oo 


X 


• 




X 


f*XJ— 2 


X 


> 




2 


<Q<X 


h- 


X 




0 


Q XQ 


< 


h- 




0 


f— h“ 


0 


0 






XOO oox 




X 




0 


ox < 


00 


X 




X 




< 


X 






XOf-2 


2 


Q 




0 


2x«-hX 




0 




X 


002 


< 






00 


H-X 2X 


1— 


0 




3 


<XXX 


< 


X 






H-QO 


Q 


X 




X 


00- 2 


0 


X 




CD 


X -XI-* 


22 


H- 






H-XX- 


00 


2 




X 


XOH- - 


— •— • 


CLI 




X 


UJ^-i X 


h-H- 






0 


Q - oO 


<< 


X 




22 


M00»-H 


h-H- 


X 




0 


XMX - 


00 y) 


X 




O'-! 


<MX- 




2 




2H- 


1 — —1 — X 


Xr- 






<< 


2- Xh~ 


x< 


<r 




0 


X o2X 


HX 


1- 




ooO 


2 - X 


h“ 


00 




XX 


X Q 


X — 


0 




X 


DCCvJ 


0X2 


X 




X2 


» j 


1 — irn — 


<r' 


c 


OX 




^ X X r*« 




2 


X>-4 


*—• nOI— 


22 11 


0 


0 


dh- 


LL2 


•—* O'— 4 


2 


•-H 


< 


XX -UJ 


f-X 


<X 


H- 


0 2h» 


oa:- 2 


Q<C I — — 




< 


00— *00 


< X 


<t 3 — 


M 


X 


2 


X oox 


XoOX*-H 




X 


oxx 


X2*-0 


XX^'-^ 


f- 


X 


o»—»0 X 


ax 2 


f-0< 


< 


X 


XH-3i-l 


2*-*2^-< 


2a-OK 


--X 


0 


K-<h- <C\J 


3HX- 


Cx 00 


— ‘f— 


0 


CL f~* •— — J 0 


2< X - 


20>aLn 


< 


XC/)H-<IXC» 


H-<X 


f-X 


oX 




Q <t *-r— 4 


xoof-rn 


OOX*~*XOX 


X 


xx>xx 


XX 


•— *X OO — r-HiO 


< 


xo o< •- 


<X 1 — 00 00 — 


h-2 — 


XCO 


2 


X X 2 — < nO 2 •- *— 4<iM 


<. X— * 


•-2 


0 


f-XOXM 0 oM *-X2M 




i-iX 


*— < 


<X 3 *<j»M 


-XJ— M 


<X M 


0 


h- 


2 CC X 0 2 rH m 


2,,. < - 


0 — 'X'«-' 


JO 


c 


23UJ XM 


r-4f-|— - 


1-~X 




h- 


X3XX— 


— -<C<I •• 


221 


X 


— « 


X2<X00<0- Qoo 


ax ro 




> 


f- > CO.--IX 


iQ 2 


•~*Oc/) >>^-o 


< 


X X * q: 


— XX 0 


h— — < — « f — 4 




X 


xxxcm— ' 


0—20X 


<!• - 


X 


000 


*--*h-X^wJ- 1) 


wH- X 


h- X < in H- 0 


0 


h-H < 


x<<XO 


00 Xf- — 


< 


x:^ 


XOO X02H- 


f-2H-2X 


*—*000 


20 


X 


<•— oox<x <*-*x ca — 


X20<XX 


h- II 


2 ‘-~'2xoxa:aof-ooM 


X XXQO 


0 


M2 *-'*XXX2X << •- 


h-MXXXX 


00 


h- 


If If If If 






•— 


<xoo 


0 




m 




xxovjm 


>r 




0 


0 


XH- 






t-H 


00 


00000 




0000 


0 



46 



to 

o 



to 



a: 

UJ 

H- 

O 

< 

< 

X 

o 



LU 

> 

< 



a: 

LLI 



>-H O 

o 

LU< 



Lu: 



3 UJ 



<t •' 
HLn 


QlcD 


O 


to o 


OLU 


(NJ 


Oh- 

XU- 


o> 

< 


o 


CC ^ 


ox 


H- 




j- 

< 


O 


H- < f— 4 


tOH 


O 


•<• 

rvj I — 1 


UJ< 

oo 


— 


o 


z: 


o 


r- 




r-4 0 




KO 


1 H 



<r •‘O- 

OXO^ OHV— 0 
^OOUJ xo<<< 
H-OH- 

vO'- 2 :<x 5 :to»-^< 

ujC 

u;o:^-j-JO« oo 

aiO<<*=^o^ OU-K2 

— iHHi-Htouj 



O 

r-l 



LUH- 

xz 

oo 



O 

in 



47 





if if if if if - 










•Jf 


>- if 












K if o 












h-LU if o 












oox ^<tUJ 










•K- 


2TH- if --UJ 










* 


LU if -J d: 










* 


QoO if<X 








f-H 




UJ if oOh- 








1 




UJI— if »-H 








X 




X3 if<0 












h-CL if >3: 








M 




5: if 








1 




1^0 if K •• 












Ko if q:lu 








X 




< if o:^: 












-JO if ooo 








M 




Ooo if < •^ 












Q-—I if >o_ 








V, 




a:< i( 












LU if X< 








CO 




if < •^ 








o 


4f 


2i< oifs:<r 








4- 




HHO:Xif rvJO 








M 




Oh- if *'UJ 


— . 






rvj 


* 


ooQ-i^ --s: 








if — 




h-o:uji^ oo 


»— • 






— iii: 




Q-OifO - 


w 






CVJw. 




2 if — 


M 


■ — ' 




O-J 




OLU^lifwO 


1 


r— 1 




if < 




HHX<if oo ^ 




I 




OO 




f-Kh-i<X^ -- 


rH 


t— J 








< xifo:-- M 


+ 


w 




+ < 




O oLUif -X 1 


►-H 


UJ 




>tM> 




c?0^if--h- — « 


w- 


QC 


F-l 


OM 




LUh-UjifOCL CNJ 


M 


o 


O 1 


4 if II 








H* 




rvIrH 


■3f 


»— 1 i Jt •»'- w> Cj Nj 




w'^ 




r^JOX 


if 


•~«o^if<r •* 




ii 




if* if O 


■K- 


cdO-K-Hif h- ^ V. 




— 


O 1 


^ ©►— ' 


if 


o if oo -- 


HH 


»— • 


1— — 


fOCO< 


if 


OOXif 0.-J-— 






2T 


0-'> 


if 


O/ if X >i») r-H 


< 


O 


O w ^ 




i^- 


<N^<if o:om-- 


f- 




o M ^ 


Mif — 


if 


ajif 


o 


— 


^ r-l-- U- 


M — X 


i^ 


ox if--occf- 


O 




V. 1 -- 1 




if 


LUOOrifC'U- *'0 


X 


rH 


f— <— O 




if 


OOOiff<^ •'—*•0 


cC 


1 


IX — -- 1 if 


cMw<^5: 


if 


OOLL-Jf --0:OX 


1 


»-H 


^ © 


00>-JO 


Q-if 


CC if M < q: 






'-'M 1 1 w-f OsJ 


4 X o<t »— t 


a: if 


CLX>if -I--- 1 


r-l 


LLI 


0 if 0 X ^MO 1 


MQ^h- 0< 


UJif 


a:f-oif OOLU-- 


•f 


QC 


»-l 0^ Mh- — X 1 -f CD 


rsj'^0»-H> 


Hif 


LU z:i< •• 




O 


O' O-JO M,-h — COif 


ifO o<w. 


z:)f 


h-T^OJis h-LUO'-'^ 


^ w- 


— h- 


< fO OLU 0 1 1 — MXif 0 


r-H'*-'— » > f— 


•— ‘ic 


XOOif< r'}-<trH 


SI <»-^o^ 


h- V. --XX-— om 


o o: 




^^Oi^ -JOOh--^ 




' oX 


O oX •' rvirvj + 


W |j W 11 O 


LUif 


1— Ujifo: OQ 


cnjoqks: 


00< oOsJ II W<<^w^ 1 < 


^ —1 o 


2ir 


iij<ta^if •'h“2TO n 


II o 


II O — 


X II xc II xo OXh-h-»-<if <if 




*-<if 


^OLLif XOOX— ' 


f-HHHX- 


--- oQ 


q:^m m rof-oo l —-if o 


II bciOrH II 


Hif 


if O 


1 


II «-H II ^d:Q-r-^ r-HO-OCXX oCO 


w>»— 4 


Z>if 


h-U.<if2 ^00^ — 


XO--’ 


— II 


---II O'-'tU LUXX— -CM 1 0< 


-j<_ih-x 


Oi(* 


o _j if o i*:: z: II LU 


rH II LU-- 


r-ix Xf\jxoOrHOOa:Q:oQ 


'itCOCCCcC 


o:if 


Oh-<if xzriii— q: II 




wf-MOrHK----f- + II II II II 


-^o— 'ooo 


CO if 


o:<oif s;<^r-<o 




Q. X II II II II II 


O*— ' »— < uoh- 


^if 




K-UL-- ii 


II Xlulu<oluxo II UM x<ll<<lu 


CyOif 


Oh-<if 0CQQQoXQQOHHOi«^^a:QQ>QQi-H02 0rNKCQ0LLU00 0a:>^“^>>c: 


if 


oo>if ^ 










if 


if 




o 


o o o 


O 


if 


if 




r-4 


fo sT tn 


<\l 


if 


if if if if if 




«— 1 


r-H rH f— i 


t— i 


oouuoo 











48 









■H- 


•» 


3* 


2 












00 




■Hr 




0 


2 












UJ 1 


LU 


* 


0 


>-m 


m3 










* 


>X Q 


X >- 




0 


OCD 


X3 










-M- 


<000 2 : 


1 — Xh- 




N^-m 




<0 












I5UJLJJ 3 


03 


* 


wm 


- 0 


0 














1 - <3 




3X 


♦.J— 


H 










*«• 


-J OXCO 


< mo 


Hr 


<x 


L'^ 


<x 












C 02 H- 


XH K-H 


H- 


00 h- 


©Q 


Xm 










* 


(X 


1— <XLL 


-H- 


•-H 


00 m 


}_ J— 












X^-iZD OLU 


OX 


■H- 


<0 


XI- 


< 







0 






LUXOX3: 


m>-m*-H 




> 2 : 


•^3 








2 




-a- 


KO-LUXO 


20 Q 




•-t— ^ 


- CL 


MX 


•» 








* 


2 <o:a:_i 


-hXI- 




h* ♦•O 




MH-»m 






2 




* 


•-hqCLLH- 


H- m3m 


if- 


CCIUCD 


mo 


-ox 


0 










o oou 


30> 


Hr 




CDO 


- 21 - 


2 




00 




* 


LL 2LUX 


XO 


H- 


W)0^ 




• • m 


•— 1 




m 




* 


0 <uja,i— 


OmOH- 


H- 


< •'OO 


Ooo 


X32 


2 




3 






f— 1 /) 


XX< 


-H- 


>Xhh 


1 — 


(XO-H 


V, 




0 






ooo o 


mLL»-3 


■H- 


»«2 < 




<XX 


00 




2 




*«■ 


LUZU ->2 


CD 2m 


-H- 


x<> 


a> 


XI- 


LU 




0 






^►-«oo< 


<<*-• q: 


if- 


< ♦• ♦• 


mo 


OX*-H 


3 










ooo 2 


»—3q: 


H- 


2 < — 


1-2 


m !2 


0 




m 






2300 LU> 


2 <xm 0 


■H- 


MOO 


3m 


xm 


2 










LU X30 


moo xo 


Hr 


♦»mcx) 


X3 


X002 


0 




in 




■Hr 


30H-02 


0^00 h-m 


■Hr 




20 


Q03 






0 






C?2<UJLJLJ 


2 


■H- 


OOrH 


Om 


-«X2 


m 




00 




■H- 


u j < (X q: 3 


(X> LL>-» 




0 ^00 


ox 


l/^l— *-H 0 






X 




■H- 


(X LLCf 


m mo< 


H- 




X 


2 X2 


Ln 




•* 






LL^> UJ 


comcD h- 


H- 




00 


02<3 


0 




m 




■w* 


x<luq: 


X<2CQ 


■a- 


060 


»-H< 


00 X 2 


00 










LL3CXXU 


mh-h-OO 




Xvj "0 


3 


2 •- 


X 




00 




•K* 


03 K 


X 


H- 




>-< 


XO<2 












OULJ < 


K 2 Xh-m 


«■ 




ooo 


CD »-H 


m 










LUOXJ— 3 


<o<x 


it 


^H-OO 


2 -« 


XI-2 






•i 




■H- 


O H-3< 


oox 0 


it 


OX •» 


m< 


Om< 


t >0 










Xo: Ooo 


'jj,u- 


-X. 


r<y L» « 


3> 


xo < 










*Jr 


<LUlLX— * 


2 mo^ 




'-^Cjc 


0 






< 


% 






Q:h-oo< 


►~«(XX 22 


it 


< **0 


m 2 


00 oO 2 V ) 0 


f^{ 


CD 


»» 






< 3> 


2 X hh 


it 


1 — — »r-(>- 


X3 


m 20 < 0 < 




m 


• 




■Hr 


m2: LUO 


Xh-m 00 


it 


OOOwO 


XX 


•-H 02 «-JCD 




2 






* 


X 00 : 2 : C 


m<X 2 x 


H- 


0 1 — t UJ 2 


►-H 


oox>-Mm 


0 


0 


0 






H- LU ►-« X 3 >- »— m h- < 0 


it 


X'-'XX 


3X — 


2 30 — 2 


0 




0 




•M- 


XLL»—: 5 :omo^ 0 


it 


QTOHH3 


2 <<M 


Xm02C0 


2 


+ 


2 






oOHO *-h 200 ooOX 


■H 


••mmO 0 


►-hhhXM 


33mXCD 










-H- 


m (XQXm 


x< 


it 


•‘UJC. 


2 h- 


0 xxm + 


2 0 




2 0 






1 — 2Xm<30ooH- 3 


it 


OX-'QCOvJ 


OOXX ^ 


XKX02 


00 


p— 1 


'--00 






0^-« O^C< 00 »-«(X oo< 


it 


CO •'C'XM ocmx 0 


x< mo 0 


2vO 


1 


*-20 






CLXX< lu 3 mms: 


it 


^(Xo MO 


>2I~ 2 


XX<X 0 


— mM 


HH 


-—mM 






S-KO CL m X <>■ 0 >cx 


it 


M<nJ-3 — 0 


«• HH 


|— 3 X 1 p— (i— (3''^ 




03 






O'~»300XLL 


0 o:o 


it 




< -2 


X < >S. 


0 ^ r-^ 


it 


mo-H 




if- 


03:0 h- 


srz-JOiZ 


-H- 


0 00 O^-H < K CD LU X M 


X 0 00 32 X om 


mXuj'-^ 


>- 


•Tr 


oC> 3<m<uo 


it 


»• •'*— < — 1 X Ll-i X 0 00 


H* 00-H < CD m X 0 


oxxo 


< 


7f 


300<Q:3H UJ 


it 


KXH OCOZCJ 


2H“,-~(m 


- <-hX2 


xxm 


2 


-xm 


OC 




<00<2^~**-h 


CD02QX 


7\- 


< •.<m^(^}CJ- •'O 


-0>l— 00<x*- X 


<KH«» QC 






m i—QUjmxi— 


it 


30o:22< 


-MO 


— Xm x<x 


•»xo 


X 


-X 


LU 




xoo>Qa:a:cx5:< 


it 


X *- 


^-XM>-i 


rHOf-mm> 


«^X i*H 




— X 


2 


^r 


mujH-mXmcLLLm X 


it 


m •»}— X 00 


ITvO -0 


-r-H<X2'-' II njo II » 


• II <t<0 II 


hH 




2 0 3 X 


or>-o 


it 


2200X II 


*•1— (• 


vO'-'UJH 


♦th r\j 


— p-H 


h- 




HH>- 3mO^-i 


mo oo^-» 


it 


n^x II 


>0 0 • ' 


-- X II — 


O'—' n 


— 


vO — 


Z) 


■Hr 


h-<cxq:cq 


XX22X 


Hr 


32 •'OOK 2 


'-H-2 - 


H-02X rH 


1— 1 


IHH 


2* *— 


o 


* 


32 < UlH-H-HHm2 


H 


Oi^2 <*-^uj<»-^trkm< coxxw 


m<'^ 


w 


m< w 


q: 




0 LU (XX 


0 3 


H 


322XOOCD2H22 ol-2(/^ h-OOI-2Or0ai-2<3 


CO 




XK3I-UJK 


XXOX 


H 


<2C2*--<ujoo»-<xmoo*-h(X*-hq xx»-«xm 


m>-«c£x 


ZD 


•n- 


(£3<0<0- 


30XH- 


it 


m03»-' 2<XO00LL(XO m>-<(XXOXO0£:XOX 


00 


■H- 


3X-^XCLoO 


X<X^“‘ 


H 


xooDO<o>3:xm 


2:xmQ(nxx:^xxox2xx 




■n- 


oOH-2h-3-H 


oomx:2 


•jf 


H 0 


•« H 


CD2 










Hr 






H 


m 


in 


r-423Q 


(\J 




>tfn 










H 


2 




<02 












•M* *Jr -K- -{f v'* 


•){■ ■« «■ -if- * ■«■ -H 


0 




xx< 








OOOOOOOUOUOOOOOOOOU 


' 0 




OCJO 









49 



to 

CC 

UJ 



0 

o 



< 

to 

I— I 

<1^ 



< 

o 

o 



X 

C\J 



00 

z 

UJ 

o 

o 

<3 

UJ 

cL 



< 

to 



< 

> 



o 

O 

O 



o 

X 

cC 



UJ--- 

HH< 

-JO 



X 
UJ| — 
:zcL 

-iiij 



id in 



O-J 

o< 

• to 

h-HH 

-J <0 
O 0>0 
O'— 

f-4-JOr 

II <00 II 

xoo II 

— h-h_jO 

Q^<s: 2 O o 

t— 4 wOOrHUJh-X njU-linUJl— ww - 

HH 7C.y-<cC ^ ll 

ULL<O»-H0f:oc*^C)'-^<i‘ 
o-»>o 











OrH 






o 




X •- 






o 




Q^tn 






-■t' 




*. f» 


o 










o 




o 




>-»LL 


rH 




X 




«— » •« 






or 




OOX 






w 




*— «fO 


H" 




X 




< •' 


a: 




— 




>^4- 


o 




— 




•>• 0 


to 




rH 






< 








HH rH 


> 


o 


o 


o 


waj 


X 


o 


X 


o 


X •- 




rH 


QC 


rH 


KX 




II 


1 


II 


Q. rH 


wa:rH:s^ 




rH-3 


UJ 


ooO ll 




H- 1 


— II 


OrH^2 



H-J0^'<0 wCL-0 — O 

O • O-JO OlXOLU — 

•.^CD>4 ' H~X^O<nZI— ^ 

II ,-i^cg II rHO O otO'— o O II *-X II ll 

II II II II XH-h- II 



wCVJOX 
rH UJ -Jl- »-H< 

•»ar< •"—! 

r-^LL.^(^JCQ 
UJ ll on ll 



o 

r- 

LL 

•* 

X 

-t 

00 







s: 


C\J 










0 


-it 




h- 






00 


0 




1— t 


C 




0 


nO 






X 




0 


< 




cc 


o 




0 


UJ 




X 


X 






o 










h- 


1 




0 


xo 




'XI 


(MOO 




Ll* 


W 




LL 






Mi 


H- — 


M 




O 


X 




Oto 


UU 




O 1 


0 


5 : 








— 2 : 


— 


X 




r— 


1— 



/ ■ ' v-' ^ 

> — j 3 :ll O -J o 



ll ►-'O 

X II U- • I— 00 “0 © 

hh- 5-5'— LL tn< 

UJCL wwUJh-SlLUUJl-rHS: II OLU u 

t-iji' " llll^SoSSoO^OO^xVi luoll<o< ll 

_j!^“-3HHHH-j:sLL-j-jootoooo_jto-^oo^-^o::OQCXUJtooo 



w< w 

«100 O 

< 1-1 *-< 

oo< H 

co*— 

CO<-^ QLCC 

vfvr> 0 '-"K 

^ o ll O 
Xf-Hh- ll 
< II O-'O'-' 
oXfo:>: 

M O'-^ 

II j00OOOrH>-00 

fsj W ' » < l~ > I + CC II 



oc:uo 
0 + 00 
oo:x 

ll X II 






^o 



OsJ 



fO 



in 



in 



o o o o 

r-4 C\1 <t 



a. 

LU 

o 

2I- 

X 

OOX 

O r-H 

•* 

xO- 

q: 

ll vOw 

z:uj< 

izujh-s: 

II o*-‘oc 

Xa^ctO 

a:a.3:iL^ 

o 
n- 



50 




I 



o 



UJ 

a 

UJ 

(X 

u. 



< 

(/) 

< 

> 

UJ 

X 



2 Z 

< 

X 



X 

UJ 







— «s 


1 - 






1— 1 


< 






UJ< 


UJ 




UJ 


2:0 


a: 




z: 


^0 


0 






-J rH 








^ ♦< 


00 




♦» 




»— 4 




*— 


X r* ' 






inm:^ 


Sd •! 








W'rH 






W' 


rH 0 


in 




OOf-^ 


i/Dr- 


0-*'* 




h- h“ 00 


•'U- 


00- 


0 


r* 


UJ 


U. 0 


0 


000 


►X 


H- r.UJ 


C\J 0 


00 •' 


--og 


2:- X 


0 


c. — 




UJXhh 


00 


0 




OOI- 




^ 1 ^ 


x< 


cd 2 : 


■ 5 f 




1 — ■ 


Q-UJUJ 


0 -- 


II WW.H- 


O-X 


»-xx 


05 : 


_I_JCL 


ujin 


OOh~ 


in ID 


rH"D H- <<UJ 


0 


o>-*UJ 


^00 


— II I— 1 — i—Oa^ooooQ 


(\JX 


OOOQdU- 


0 -- o\ 


— -0 — 1 — OOOQ<»~iH-^ 


— oC 


vO U-O 


C*if\xu^ 


oooQO »-<<— 


— J 


•Jc 




OQ Ilo 0>>0 


OU- CD 


<0 


♦'iic^rooOw 


OOi— IrH II II II 0 oO' 


1—4 »■ 


siin- • 


r— 1 UJ 00 rH 


O'*-" 0 0 II "^11 H—lXJr-^ 


♦'X II 


0 »• •* 


II + II o^COXHI— '-'O O-J ».invOr-i 


UJH- 4 s 0 » l-Hl 


X c\j 


»-H_jO''~'in 0 ovO>t 




Z)00-"W 0: 



r-H— 3 II II u_ o orHCNjini'-.-i-^oo-- H -t X z: 

O Xlui O ►-h->-^w.ww^w^w00c/)UJOUj<whh II ulkhjLhhO: 

Ch II w.w.UJUJLULLlUJUJw.wH-h-»--:?:UJK h-x -i~x 

II hh :zhq 

Oi^rHO-wo^ u- u.*-»H-4»-Hi-Hi-H»-HULa.a:oa:o»-^o>-^a'o oujz: 

Q :^ (/) HHC/) O CO HH HH -j -j _j i-i (-H 3 : o 3: LL -J o 00 3 : LL > o cr UJ 

O LO 



00 

csjcvj 



in inomc- 
ro 



in 



51 



= 




SUBROUTINE RUNK 



•• 



•K 

4f 

* 

rr 

•K- 

* 



I 00 

HH Qh-4 

-J UJ 

•-•ctlh-^ 

KLU<0 

3QH-.H-. 

a:oh« 
QOO< 
2 003 

<Q00O 
2<LU 
OO 

aio>*Lu 
^-UJ^-X 
C OO^-*!— 
UJ o 
o:ajOLL 
OX-JO 
HUJ 



•if 


> >5: 


2 


•if 


-joo ct: 


2 




OOLU^JU 


3 


•If 


3XX 


— J 0 


V 

A 


O-JO 


QX 


•if 


H-tO*— 'LU 


o^- 


4f 


>OOf-X 


a 


•if 


LU q:i— 


o:lu 


•if 


a:oLu 


LUQ 


•if 


CLO> 0 ^ 





o:3lu 
^OOfO 
2 UJ< 

ox 00 
ct:i— — Joo 
g:<< 
aj3*-^Q- 
20H- 
i-iLLz:uj 
K LUX 
3<q:i- 
o LU 
q:ollx 
cD^ru-i- 

0»-H 

O0MQ3: 



•if ->f -if -if -if 


rf 




LUX 




•if 


— <r 


CD3 




•if 


O OsJ 


-J 


O 


•if 


o -- 


H- <t 


1-2 


* 


V|-LU r-H 


oo> 


LU 


•if 


'-LUX 


3 


t/)0 


* 


-JXCL 


5IX 


1— H-4 


•if 


<x< 


> 


O.LU 


•if 


001— a: 


•— M— 4 


2 


•if 


^o 


XK 


UJ2 




<o •* 


CL< 


1-LU 


•if 




O 


f-K 


•if 


-h-r- 


XX 


< 


•if 


f— •'vO 


02 


h- 


■«• 


q:x'— 




200 


•if 


02UJ 


x< 


<q: 


•if 


ooox 


00 


q:^ 


•if 


< *>DC 


200 


ox 0 


•if 


>CLX 


XX 


O 2 


If 


-2k 


002 


XXO 


■if 


X <3: 


3 


CLXh-c 


•if 




XOO 



LU < 
Xf-Z) 
l“20f 
LULU 
0000 
LUlU-J 
^oc< 

OCLH-. 



“«• 

* 



Mor^ 

•*LUs0 

ooo 

O ••5: 

--03: 
OO 



XOO 

f-< 

20 

occ: 

I— i<i 

1-3: 

<2 

3^ 



•if 


X LU 


•if 


^1—0 Q 0 


2X1- 


•if 


XM--- OO 






•if 


00KX> 




xQ 


xqc:2 


•if 


o:--h- XQ 






•if 


xxx< 


CM 


:sx ocx:x 


3 X 


■>f 


-XsO 






•if 


-J2 1— 


OO 


X X- >- < X 


ox a: 


if 


— 1 oox 






'A 


jO 


i 


X ^ 2 O / '^ 


' 1-* *3 ‘ ' - 


-X 


OQ-UJ--‘- (jO 






•if 


<<ox 


CM 


KK3203 


DCt-iX 




r» ^ 






•K- 


-Il-2<1 


2 


Hi^XooO 


XXX 


if 


— ooo<cx 






•if 


q:|---«2 




200 3 OO 


3: •— 1 


if 


< ^3z<t\-cx: 






•if 


<32ct: 


•if 


— i^XQf>- 


Q Q 


•if 


1——' ►wCOO 






•if 


>^ct:x 


— - 


X>XO> 


XOO 




00 0'^'“‘02 






•if 


1 luk 




>>Q<ct:20 


klux 


if 


O »— 4 • — 1 X 1 H- 1 




o 


•if 


XX>2 


CM 


K ISXX2 


oox 


if 


X-^rHXX 




f-H 


•if 


XOO-1 


X 


^U- 3X 


XOK 


•K- 


-x> 






•if 


H20 


1 


OOXXCr3 


x2 


■>f 


••XOO—' 1— 




o 


•i;- 


3 2 


CM 


o xxxc:r 


X X 


•if 


^Qc: 00 O' 2^-4 




1- 


•if 


ooo:2< 


00 


xxl— Kct: X 


002K 


if 


O X X nJ* 00 






■if 


X o 




x> u Ct: 


Xh- 1 


•if 


CO -X)-- 2 




Q 


if 


ildrar'-iLL 




>H-4U_X X 


2K3: 




mxx 




O 


•if 


<ujkq 


CM 


KOOX 


X 


Vv* 


M< •'X3Q 


— 




■if 


1-Q< 




x< << 


kkq 


•if 


••K — ' X 


•— 1 





<>LULUk-._J 

o^ci:d>-»< 
HH Ol<KQT OO 
1— lu3 3lu^— 



00 LU 
UJQO- 
X*-i< 



m 




o:qoo2< 


o 


y: 




X o0O0»--4> 


XXO 


•if 




>X 


OXOO 


>— 4 




xxxxx 


1—00 


X 




XKXXXX 


X < 


X 




X Hl-l-1- 


ox 


1 


t • 


1—00 


<1-00 




X 


•— loooooooo 


X3X 


II 


O' 


CO »— IHH l-H 


X3 




Ul 


HH *-H 


X2X 


z 


X 


XCMCMCMCM 


oo< 


y: 


3 


yiCL^LOlLZ 


xo> 



-n- 

k 



4f 
■«• 

•if 
•if 
•if 
•if 

•if ^ -if -if -if •«••«• -if -if •«• -if if -if ^ -if -if -^f -if -if -if -if -if ^ -if 



O 

•—I — ‘i— - 

OoOr-lctoOLL o oOO 

•• •<?ct:^-»oo-J»-H^-‘oo 

I— LU(\lO 0< ovO-^ 

<x •'OOLU'4'OOfHOo*;;- o 

— JOI2--'^l— LULU<UA 

Q: ► 0'X<r-H< II cc: OSIr-i 

•«i-2vt<Q: II > u.<o - 

SOO'-^ I »-H I I— I ^HHsO 

OH-cf\JM< II II OOO'^ 
2 •'OOX *-H 

OXZICLUJI— rHi-tC<00 II LU 
22LU<X<CvJ'--'r-i5:-^ I— 

2<2a:f~x »-H o 

o-J^o f-oxo'-<x»-Ha: 

OCDQ **LU Q0.O00HH002 

•if 



OUOOOOOOOOOOOOOOOOOOOOOOOOO 



if ox 




20 


O' 


^ZD 


r—i 


0000 


CM 


OO 





52 



150 FORMAT ( ’ 1 • , 20 X , • C OMPUT AT I ON OF FIRST THREE MODES FOR FREQUENCY 
='!'F8o5,' CYCLES/MIN, • t/ ) 
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DO 301 M=l,199 

3C1 IF (VHM)=:^W(M+1 ) oLToOo ) KK=KK+:i 

IF ( ABS( ENOVAL )o LToEPS) GO TO 350 
315 IF(LL<,EQ«1) GO TO 330 

3C6 IF (ABS(W(201) )oGToABS(W(2C0) ) ) GO TO 308 
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